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ABSTRACT 

It is shown using Lund Monte Carlo that, unlike the average prop- 
erties of the hadronic system inside jets, the anisotropy of dynamical 
fluctuations in these systems changes abruptly with the variation of 
the cut parameter y cut . A transition point exists, where the dynamical 
fluctuations in the hadronic system inside jet behave like those in soft 
hadronic collisions. Thus the anisotropy property of the dynamical 
fluctuations can serve as a probe for the soft and hard processes in 
high energy collisions. 
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As is well known, the presently most promissing theory of strong interaction — Quan- 
tum Chromo-Dynamics (QCD) has the special property of both asymptotic freedom and 
colour confinement. For this reason, in any process, even though the energy scale, Q 2 , is 
large enough to be able to do perturbative QCD (pQCD) calculation, there must be a non- 
perturbative hadronization phase before the final state particles can be observed. Therefore, 
the transition or interplay between hard and soft processes is a very important problem. 

In current literature, this transition is determined by some cut-parameter. For example, 
in doing theoretical calculation a parameter Ql is introduced. When Q 2 > Ql the pertur- 
bative QCD is assumed to be applicable and the process is hard. While when Q 2 < Qq the 
perturbative calculation is unallowed and the process becomes soft (nonperturbative). How- 
ever, the value of Ql is not determined exactly. It decreases steadily as the developement of 
perturbative technique. 

In experimental data analysis people use some "jet-algorithm" (e.g. Jade or Durham |J 
ones) to combine the final-state particles into "jets". Each jet is assumed to be originated 
from a hard parton, and the hadrons in the jet is produced softly from this hard parton. 
Thus the transition between hard and soft processes is described as the production of hard 
partons and the subsequent hadronization of these partons. In this formalism there is also 
a parameter — y cnt . The value of this parameter determines how the hadrons are grouped 
into jets, and whether an event is a "2-jet event" or a "3-jet" , "4-jet" ones, 

Let us concentrate on the 2-jet events. By definition, these two jets should be developed 
softly from two hard partons and no hard process is involved in the evolution. If there is 
any hard process in the developement then we say that a third jet appears. Historically, 
it was the observation of the third jet in e + e~ collisions that confirmed the existence of 
gluon |§-||. In this sense, there should be a definite value of y cut , which is consistent with 
the physical meaning of "jet" . 

On the other hand, due to the success in pQCD calculation of jet, people sometimes 
take the number of jets in an event as indefinite, depending on the value of y cu t, which can 
be chosen arbitrarily. Their stress is in ultilizing this dependence to confront the pQCD 
calculation with experiments. From this point of view, the physical meaning of jet and the 
associated concepts - - "soft" and "hard" are neglected. A process is hard or soft is not 
judged physically, but is determined through the technical problem of whether the process 
can be calculated by perturbative QCD. 

Let us remind that physically, soft and hard are distinguished through the magnitude of 
transverse momentum. In hadron-hadron collisions at energies below top-ISR most of the 
final-state hadrons have low transverse momenta and the process is soft. At collider energies 
high-transverse- momentum jets, coming from hard parton collisions, start to appear M. The 



to as semi-hard. The critical value of transverse momentum for the transition between soft 
and hard (semi-hard) is about 4-5 GeV/c. 

The following important questions arise: 1) Does the number of jets in an event possess 
any definite meaning? If yes, how to determine this number, i.e. how to decide the correct 
value of 7/ cu t for the determination of this number. 2) Is it in principle possible to locate the 
transition between soft and hard processes in the hadronic final states of high energy e + e~ 
collisions? If yes, how to do that? 

In order to answer these questions, let us remind that the qualitative difference between 
the typical soft process — moderate energy hadron-hadron collisions and the typical hard 
process — high energy e + e~ collisions can be observed most clearly in the property of dy- 
namical fluctuations therein. It is found recently || that inspite of the similarities in the 
average properties, the dynamical fluctuations in the hadronic systems from these two pro- 
cesses are qualitatively different - - the former is anisotropic in the longitudinal-transvere 
plane and isotropic in the transverse planes while the latter is isotropic in three dimensional 
phase space. 

This observation inspired us to think that the dynamical-fluctuation property may pro- 
vide a probe for the transition betweem soft and hard processes inside the hadronic final 
state of high energy e + e~ collisions. In the present letter we show, using Lund Monte Carlo 
simulation, that this is indeed the case. 

In total 500 000 events are generated for 91.2 GeV e + e~ collisions using JETSET7.4. 
The resulting hadronic systems are analysed using Durham and/or Jade jet-algorithms. The 
fractions R 2 , R3, -R4 of the 2-, 3-, 4-jet events in the whole sample are plotted versus the 
value of y cut for both Durham and Jade algorithms in Fig.l. It can be seen clearly form 
the figures that the definition of "jet" depends strongly on the value of y cut . When y cut 
is big, most of the events are taken to be "2-jet" events. In the limit of very large y cu t, 
the whole sample consists of only "2-jet" events. On the contrary, when the value of y cut 
decreases continuously, the jets are divided further and further, and gradually most of the 
events become "multi-jet" (more than two jets) ones. 

At the energy in consideration, it is certainly impossible that all the events are 2-jet ones. 
Neither is it possible that most of the events are multi-jet ones. In order to determine a rea- 
sonable value of y cut , we have to use the dependence of some physical property of the system 
on y cut . As example, we show in Fig. 2 the dependence of average charged multiplicity {N ch ) 
and average ellipticity (e) on y cut for the "2-jet" sample determined by Durham algorithm. 
The ellipticity e is an event-shape parameter defined as the ratio of minor T 3 to major T 2 in 
thrust analysis |]TD|H 

e = T 3 /T 2 . (1) 



expected that, when y cnt increases, more and more "inpurities" (multi-jet events) are mixed 
into the "2-jet" event sample, and the jet cone will diviate more and more from being circular. 
So the average ellipticity (e) will drcrease with the increasing of y cnt . It can be seen from 
the figure that this is indeed the case. However, the value of (e) changes smoothly with y cut 
and it is hard to get a probe for a reasonable value of y cut by using (e). The same holds also 
for (iV cll ) and other average quantities. 



r Durham 






A 


; Rj/ 













_- uade 




Fig.l The ratio of 2-, 3-, 4-jet events 
as function of y cu t 
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Fig. 2 Average charged multiplicity and ellipticity 
as function of y cut 



Let us turn now to consider the dynamical fluctuations. These fluctuations can be char- 
acterized by the anomalous scaling of factorial moments (FM) ||11|| : 
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where a region A in 1-, 2- or 3-dimensional phase space is divided into M cells, n m is the 
multiplicity in the mth cell, and (• • •) denotes vertically averaging over the event sample. 
Note that when the fluctuations exist in higher- dimensional (2-D or 3-D) space the projection 
effect |L2 will cause the second-order 1-D FM goes to saturation according to the ruleQ: 
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where a = 1,2,3 denotes the different 1-D variables. The parameter 7 a describes the rate 
of going to saturation of the FM in direction a and is the most important characteristic for 
the higher- dimensional dynamical fluctuations. If 7 a = 7& the fluctuations are isotropic in 
the a, b plane; while when 7 a ^ the fluctuations are anisotropic in this plane. The degree 



of anisotropy is characterized by the Hurst exponent H abl which can be obtained from the 



values of j a and 7^ as [14 



= ^ (4) 

The dynamical fluctuations are isotropic when H ab = 1, and anisotropic when H ab 7^ 1. 



For the 250 GeV/c 7r(K)-p collisions from NA22 the Hurst exponents are found to be [15 



H PtV = 0.99 ± 0.01, H ypt = 0.48 ± 0.06, H w = 0.47 ± 0.06, (5) 

which means that the dynamical fluctuations in this moderate energy hadron-hadron col- 
lisions are isotropic in the transverse plane and anisotropic in the longitudinal-transvere 



planes. This is what should be JT6[ , because there is almost no hard collisions at this energy 
and the direction of motion of the incident hadrons (longitudinal direction) should be pre- 
vileged. Note that the special role of longitudinal direction in these soft processes is present 
both in the magnitude of average momentum and in the dynamical fluctuations in phase 
space. 

In high energy e + e~ collisions, the longitudinal direction is chosen along the thrust axis, 
which is the direction of motion of the primary quark- ant iquark pair. Since this pair of quark 
and antiquark move back to back with very high momenta, the magnitude of average mo- 
mentum of final state hadrons is also anisotropic due to momentum conservation. However, 



the dynamical fluctuations in this case come from the QCD branching of partons |17| , which 
is isotropic in nature. Therfore, although the momentum distribution still has an elongated 
shape, the dynamical fluctuations in this case should be isotropic in 3-D phase space. 

A Monte Carlo study for e + e~ collisions at 91.2 GeV confirms this assertion ||. The 
dynamical fluctuations are approximately isotropic in the 3-D phase space, the corresponding 
Hurst exponents being 

H Pt(p = 1.18 ± 0.03, H vpt = 0.95 ± 0.02, H w = 1.11 ± 0.02. (6) 

The present available experimental data for e + e~ collisions at 91.2 GeV also show isotropic 
dynamical fluctuations in 3-D |T8| . 



Now we apply this technique to the "2-jet" sample obtained from a certain, e.g. Durham, 
jet-algorithm with some definite value of y cut . Doing the analysis for different values of y cut , 
the dependence of dynamical-fluctuation property of the "2-jet" sample on the value of y cut 
can be investigated. 

Let us try to discuss what results can be expected? 

As we have shown in Fig.l, when y cut is very big the "2-jet" sample coincides with the 
whole event sample, R2 = 1. In this case, the fluctuations are known to be isotropic in the 
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Fig. 3 The variation of 7 with i?2 and y cut 

As the decreasing of y cut the multi-jet events, which contaminate the "2-jet" sample, 
will be cleared away gradually, and at a certain value of y cut , a "pure" 2-jet sample will be 
formed. The word "pure" is used here to indicate that these two jets are developed softly 
from initial partons and no other jet(s) has been mixed in. 

It can be expected that the dynamical fluctuations in the "pure" 2-jet sample will mimic 
those in the soft hadronic collisions, i.e. isotropic in the transverse plane and anisotropic in 
the longitudinal-transverse planes (7^ = 7 V 7^ 7 y ). 

Thus the variation of 7's with the decreasing of y cut (or decreasing of R 2 ) is expected to 
be: At first, when y cut is very big, the "2-jet" sample is identical to the whole event sample 
(R 2 = 1), and the three 7's equal to each other; As the decreasing of y cnt (the deceasing 
of R 2 ) the three 7's depart, and becomes, at a certain value of y cut , isotropic in {pt,f) and 
anisotropic in (y,p t ) and (y,tp), % t = 7^ 7^ 7 r 

The results of simulation are presented in Fig. 3 (a). It can be seen from the figure that 
the above expectation comes true. The characteristic behaviour j Pt = 7^ 7^ 7 y arrives at 
1/cut ~ 0.0048 (R 2 ~ 0.48). The values of 7's and the corresponding Hurst exponents at this 
point are listed in Table I. For convenience we will call this point, where 7 Pt = 7^ 7^ j y , as 
transition point. 



Table I Parameter 7 and Hurst exponents at the transition point 
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Fig. 4 Comparison of the speed of going to saturation of 
F 2 for different 1-D variables at different R 2 



Note that in the Durham algorithm that we are using the test variable y is essentially 
the relative transverse momentum k± squared [[IS]. The transition point y cnt « 0.0048 
corresponds to ~ 4 GeV/c, which is consistent with the critical value of transverse 
momentum between soft and hard (semi-hard) components in hadron-hadron collisions. 

It is instructive also to follow the evolusion of 7's with the increasing of y cu t (incresing 
of R 2 ). 

It can be seen from Fig.3(a) that, when y cut (R 2 ) is very small, where the two "jets" are 
highly undeveloped and each consists mainly of one hard parton, 7^ is consistent to zero, i.e. 
there is no dynamical fluctuation in (p at all. On the other hand, at this point 7 Pt is almost 
as large as j y , showing that the dynamical fluctuations in this undeveloped "2-jet" system 
behaves as an isotropic 2-D fractal in the (y,Pt) plane. 

When y cut (R 2 ) increases, j Pt departs with 7 y and approaches to 7^. What is important 
is that 7 Pt and 7^, instead of going up parallelly, cross over each other, turns from 7^ < j Pt to 
lip > 7 Pt 5 resulting in a sharp transition point. After that, the three 7's approach eventually 
to a common value, and the "2-jet" sample approachs to the whole event sample. 

In order to show the evolusion of the anisotropy property of dynamical fluctuations with 
the variation of y cut (R 2 ) more clearly, we take three typical points: (A) R 2 = 0.18, (B) 
R2 = 0.48, (C) R 2 = 1, indicated by arrows in Fig.3(a). Point A corresponds to the 
case of undeveloped jets, B is the transition point and C is the whole sample. Since the 
anisotropy property of dynamical fluctuations determines solely by the rate of approaching 
to saturation of FM, which is characterized by the parameter 7, we rescale the F 2 (pt) and 
F 2 ((p) appropriately, letting them coincide with F 2 (y) at M = 3 and arrive at a common 
saturation height with F 2 (y). The results are shown in Fig. 4. It can be seen from the figure 
that when R? = 0.18, F 2 (lp) does not increase with M, i.e. no dynamical fluctuation at all in 



R 2 = 0.48 (transition point), F 2 (p t ) and F 2 (ip) go to saturation almost with the same speed, 
much slower than F 2 (y) do. When R 2 = 1 (whole sample) all three F 2 coincide and go to 
saturation with an identical speed. 

For comparison, we have also done the same analysis using Jade algorithm. The results, 
shown in Fig.3(6), are qualitatively the same: At small y cut (R 2 ), 7^ vanishes and j Pt ~ 7 y ; 
As y cut (R 2 ) increases 7 Pt and 7^ approaches each other and cross over at y cut ~ 0.158 
(R 2 ~ 0.39). This is the transition point for Jade algorithm. The parameter 7's at this 
point are j y = 1.22 ± 0.04, 7pt = 0.51 ± 0.09, 7^ = 0.59 ± 0.08. 

In this letter we have shown using Lund Monte Carlo that, unlike the smooth change 
of average properties of the hadronic system inside jets, the anisotropy of dynamical fluctu- 
ations in these systems changes abruptly with the variation of the cut parameter ?/ cut . At 
yfs = 91.2 GeV, the dynamical fluctuations in the whole e + e~ collision sample (large y cut 
limit) are fully isotropic in the 3-D phase space, and become highly anisotropic (almost no 
fluctuation at all in ip) for small y cnt where the "jet" is highly undeveloped. A transition 
point exists, where the hadronic system inside jet behaves like that of the soft hadronic col- 
lisions, i.e. the dynamical fluctuations are isotropic in the transverse plane and anisotropic 
in the longitudinal-transverse planes. The corresponding relative transverse momentum at 
the transition point is about k± ~ 4 GeV/c, which is consistent with the critical value of 
transverse momentum between soft and hard (semi-hard) components in hadron-hadron col- 
lisions. Thus the the transition point determines the physically meaningful value of y cut , and 
thereby gives the number of jets in an events. The anisotropy property of the dynamical 
fluctuations can serve as a sensible probe for hard and soft processes. 

This observation is not only meaningful in the study of jets in e + e~ collisions but also 
enlightening in the jet-physics in relativistic heavy ion collisions, which will become impor- 
tant |20| after the operation of the new generation of colliders at BNL (RHIC) and CERN 
(LHC). 
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